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^ . Abstract 



We study the dynamics of Einstein's equations in Ashtekar's variables 
from the point of view of the theory of hyperbolic systems of evolution equa- 
tions. We extend previous results and show that by a suitable modification 
of the Hamiltonian vector flow outside the sub-manifold of real and con- 
strained solutions, a symmetric hyperbolic system is obtained for any fixed 
choice of lapse-shift pair, without assuming the solution to be a priori real. 
We notice that the evolution system is block diagonal in the pair (cr a , At), 
and provide explicit and very simple formulae for the eigenvector-eigenvalue 
pairs in terms of an orthonormal tetrad with one of its components pointing 
along the propagation direction. We also analyze the constraint equations 
and find that when viewed as functions of the extended phase space they 
form a symmetric hyperbolic system on their own. We also provide simple 
formulae for its eigenvectors-eigenvalues pairs. 
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1 Introduction 

A substantial effort has been undertaken in the last decade to understand 
what a canonical quantization of space-time would mean in terms of self- 
dual variables and the corresponding loop representation it implies. Most 
of the effort has been geared to the study of the constraint equations and 
their algebra, while little has been done about the evolution equations, in 
the understanding that, since the Hamiltonian of the theory is a linear com- 
bination of the constraints, the information contained in their algebra is all 
what is needed. 

Nevertheless we believe that the study of the evolution equations as such 
should be of relevance for a better understanding of the constraint algebra, 
for it might assert that it has very important properties which might other- 
wise be very difficult to recognize, namely that of giving rise to a well posed 
initial value formulation, that is, a theory in which one could predict the 
future based on knowledge gathered at an initial Cauchy surface. 

It is not clear a priori that a classical well posed initial value formulation 
is essential for quantization. However all physically interesting quantum 
theories we know of correspond to well posed classical systems, and this fact 
is implicitly used in most approximation schemes. 

The well posedness of the initial value problem for the classical theory in 
the usual tensorial variables was asserted in the fifties by Choquet-Bruhat Q 
using a particular reduction of the equations in the harmonic gauge. Since 
then, a number of other formulations of the theory as symmetric hyperbolic 
systems have appeared (see H, |, |, |, Ill, Q and for an updated review 
with a more complete reference list see |8| ), mainly aiming to gain some 
advantages in dealing with certain specific problems, such as asymptotic 
studies, global existence, the Newtonian limit and numerical simulations of 
fully general relativistic configurations. 

It is not clear that the study of the dynamics of Einstein's equations 
in Ashtekar's variables as a symmetric hyperbolic system would have an 
important impact on applications such as those mentioned above, but its 
intrinsic beauty, simplicity, and economy can not be in vane. 

In dealing with first order quasi- linear systems of evolution equations, 
-as Einstein's equations naturally are in Ashtekar's variables- a sufficient 
condition for well posedness is symmetric hyperbolicity, that is, if the system 
can be written in the form: 




(1) 
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where u is a "vector" and M°, M a and B are matrices which have smooth 
dependence on u and, in addition M° and M a are symmetric and M° is 
positive definite []. This condition is the one of more immediate use, in 
particular for numerical simulations, for it enables one to estimate the growth 
of the solution in terms of conventional norms; therefore it is always worth 
checking whether or not it holds. If a system does not have this property it 
does not mean that its associated initial value formulation is not well posed, 
for the necessary condition for well posedness is a weaker condition, namely 
strong hyperbolicity 0, that is, that the eigenvalues of the principal symbol, 
iM a (u)k a be purely imaginary and that it has a complete set of eigenvectors, 
all of them smooth in their dependence on u and k a . It turns out that the 
Ashtekar's system can be extended outside the constraint and reality sub- 
manifold of the phase space in such a way that both hyperbolicity conditions 
hold. 

In determining whether a formulation of general relativity is well posed or 
not there usually arises the problem of constraints, since we can not locally 
separate the fields into free and into dependent ones with respect to the 
constraints. We have to think on the evolution equations not only as giving 
the dynamics at the constraint sub-manifolds, but also as evolution equations 
for the complete fields, that is, as equations also valid in a neighborhood of 
the constraint sub-manifold. But this makes those equations non-uniquely 
determined; all dynamics which at the constraint sub-manifold give the same 
evolution vector field are valid, and therefore equivalent. Thus, if we modify 
the equations, by adding to them terms proportional to the constraints, 
we obtain equivalent systems of equations. Thus, in order to show well 
posedness one only has to search for one of these equivalent systems of 
equations in which the above conditions hold. 

In the Ashtekar formulation of general relativity a further complication 
arises because besides the constraints there are reality conditions. Along 
evolution the soldering forms are in general not anti-hermitian, even if they 
started that way at the initial surface, so without assuming very stringent 
evolution gauge conditions or imposing elliptic gauges, one can only require 
the metric reality conditions along evolution and since these conditions are 
non- linear in terms of the basic dynamical variables, it becomes necessary 
to deal with them as if they were a new set of constraints. 

In a recent paper ||, the authors have answered the question of 



1 For a short introduction into the topic see Q and references therein. 

2 The name strong comes because it is stronger than the primitive concept of hyperbol- 
icity, namely that the eigenvalues of the associated symbol be purely imaginary, but it is 
a weaker condition than symmetric hyperbolicity. 
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hyperbolicity by finding a representative within the equivalent class of evo- 
lution equations and showing that it is a symmetric hyperbolic system. In 
that paper it was assumed that the soldering form was anti-hermitian, re- 
stricting unnecessarily the set of gauge conditions (choices of lapse-shift pair 
and of the SU(2) rotation gauge) for which symmetric hyperbolicity holds. 
That restriction can be easily lifted to allow for all soldering forms to be 
SL(2,C) rotated to an anti-hermitian one, that is, for all soldering forms 
satisfying the metric reality condition. We do this in the next section, (§^). 

Unfortunately this lifting is not enough to grant the well posedness of 
the system, because a small departure from the reality sub-manifold would 
imply an immediate and strong instability for certain perturbations. Thus, 
the usual methods of constructing the solution via a contraction map of 
successive approximations can not be used, nor can one use these equations 
for numerical simulations. To overcome this problem a further modification 
to the evolution equations is needed, this time outside the reality condi- 
tion sub-manifold, which makes the system symmetric-hyperbolic in a whole 
neighborhood of that sub-manifold. We study this modification in §||. In or- 
der to distinguish it from the prior modifications above mentioned, we shall 
call it Bow regular ization. 

We then study, in §||, the diagonalization problem for the principal sym- 
bol of the system, that is, the issue of strong hyperbolicity and of the prop- 
agation speeds of perturbations. In particular we display remarkably simple 
expressions for the different eigenvalues-eigenvectors of the system. We con- 
clude that the system can not be extended outside the reality condition 
sub-manifold in order to make it strongly hyperbolic there, keeping at the 
same time unchanged the simple dependence of the characteristic eigenvec- 
tors on the solution. On the other hand one can modify the equation into 
a strongly hyperbolic pseudo-differential system, the advantage here being 
that the characteristics of the system remain very simple. 

In §H| we study the issue of the constraint and reality conditions propa- 
gation, showing that the system stays on those sub-manifolds if initially so. 
We show in particular that the constraints propagate also via a symmetric 
hyperbolic system. The diagonalization of this system is also remarkably 
simple. 



2 Symmetric Hyperbolicity 
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2.1 Earlier Results 

In a recent paper (9j , the authors have shown that the Einstein's equations in 
Ashtekar's variables constitute a symmetric hyperbolic system. In that paper 
it was assumed that the soldering form was anti-hermitian 3 , restricting 
unnecessarily the set of gauge conditions for which symmetric hyperbolicity 
holds. 

In order to obtain a symmetric hyperbolic system for any given lapse- 
shift pair, we suitably extended the field equations outside the constraint 
sub-manifold. Since this extension only involves the addition of terms pro- 
portional to the constraints, we know that the evolution vector field is tan- 
gent to the constraint sub-manifold, and so the physically relevant evolution, 
that is, the dynamics inside this sub-manifold, remains unchanged. 

For completeness, we repeat here part of the calculations made in ||, 
showing in addition that the above extension is unique. 

It is known ((nj, see appendix ||| for a brief overview) that, using the 
Hamiltonian formulation to determine the dynamics of the pair (a a ,A b ), 
Einstein's equations can be written as 



Cfa b 




C t A b 


= C^A b + ^=N[a a ,F ab ], 


C(a,A) 


:= tr(d a a b F ab ) = 0, 


C a (a,A) 


:= tv(a b F ab ) = 0, 


C A B (a,A) 


:= V a a a A B = 0, 



Since the evolution equations have a block diagonal principal part, that 
is, in the time evolution for a a (A a ) there only appear space derivatives of 
a a (respectively A a ), the desired modification is of the following form: [] 



C t a b = C^a b -^=V a (N[a a ,a b })+ a [C,a b ], (3) 

3 Since we are using intensely the matrix algebra of the soldering forms we refer to them 
as anti-hermitian in the matrix sense, in the abstract sense, as denned in jlo[ they are 
actually hermitian. 

4 One can always add to these equations given SU(2, C) gauge terms in the usual way 
without affecting hyperbolicity, for if they do not depend on a" , nor on A a , they do not 
enter the principal part. 
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C t A b = C^A b + -j=N[a a ,F ab }+Pa b C + ie b dc d c C d . (4) 

where a, (5 and 7 are complex functions to be determined so that the system 
becomes symmetric-hyperbolic, meaning that the principal symbol becomes 
anti-hermitian with respect to the canonical inner product 

(1*2,1*1) = ((^ 2 ,^ 2 ),(ai,A 1 )) 

= ti{d\ a a\)q ab + tT(A^ a A\)q a \ 

and considering the soldering forms, a a anti-hermitian. 

Recall that the principal symbol of a quasi-linear evolution equation sys- 
tem 

u l = B l j a {u)V a u j +M l (u) 

is given by P(u,ik a ) = i B l j a (u)k a . In our case u denotes u = (a a ,A b ). 
Thus we need to prove that 

P12 + P ] \2 = {u 2 , Pm) + (i* 2 , P^m) = 0. 

Using the fact that any complex matrix can be written as u = — -^tr (ua e )a e 
and that [<7 a ,<r b ] = a\pie ah c a c , the principal symbol can be written as 



k N a k 
Prfaika)*!" = -i-2-5-tr(ai b £ e )a e - Ntr^ 6 )^^ 



(a - _N ) katifa a a e )e eb d a d , (5) 



k N a 2 
P A (a,ik c )A\ = -i^^ t r(A lb d e )a e + -m [ MA 1 b] a e )e ae c a c 

+i pV2ae cde k c a b tr(a e A\) + i2 je b dc a c tr(a a k [d A a p) 

A very simple calculus shows, using the inner product above and the anti- 
hermiticity of the background soldering form, that 

k N a + k N a 

PcYi + PKvi = -i^^tr(ai 6 a e )tr(^ 6 a e ) + i^-^tr(a 2 ^ b a e )tr(a lb a e ) 

-N k a tr(ai 6 a e ) £ ae c tr(4 b a c ) - -N k a tr(d lb a e )e ac e 



a a 

n/2 



xtr(a 2 t fe a c ) - !X_( Q - 4=N) fc a tr^e)^ tr(4 b a d ) 



(a + -=N) A; a tr(a 16 a a ) £ e d tr(,r 2 Ta a e ). (7) 
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The sum of the first and the second term vanishes because of the symmetry 
of the background metric used in the inner product, the sum of the third 
and the fourth term vanishes because of the anti-symmetry of e abc , finally 
the last two terms should vanish in order to get the desired result, but it is 
impossible for the sum to vanish because there are different contractions on 
the indices and o\ is independent of 02, thus, the only way for these terms to 
vanish is that each one of them do so separately This is obtained requiring 

a = -4=N. 

Next we calculate 



h AT a h N a 

Pau + P f Ai2 = -i^tr(AV e )tr(A t2b cr e ) + ,^-tr(A 16 a e ) 

xtr^Ve) +N — £ aem tr(AVe)tr(At 2b a m ) +N — e aem tx(A lb a m ) 
u a 

tr(AV b a e ) + i(3V2ae cde ti{^ 2b d b )k c tr(A Ve) - ijk d e dc b tv(A^ 2a a a ) 
k c ti(A lb d c ) + ije dc b k d tr(A^ 2b a c ) ti{A la a a ) - 0V2ae cde k c tr(A^ d a e ) 

tr{A lb a b ) -N ^e apn ^ 2b ~a n ) tr(A 1 a a p ) - i^e a dp k b tv(A^ 2a d p ) 
tT(A 1 d a p ) - N — £ a P n tr(AV a a p ) tr{A lb d n ) + i-fe b dc k a tr{A^ d d a ) 



a 



tr(A lb a c ). (8) 

The sum of the first and the second term vanishes because of the symmetry of 
the metric, the third and fourth terms vanish because of the anti-symmetry 
of e abc . The next four terms give rise to 

i tM 2]b °b) tr(A VeK* + 1) + i tr(A 2 V e ) tx{A l a d a )e dce k c 

x (pV2a + 7) 

and there are no other terms of the same kind, so they must vanish. The 
only way (since A 1 and A 2 are independent) this can happens is if 

(3y/2(T = -7. 

Finally the last four terms vanish if we use 

2A [ a 1 = £ EdmeA a , 

and choose 

% 1 

7 = N, and consequently (3 = ;r~/=N- 

c (j 2 v2 
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2.2 Lifting the Anti-Hermiticity Condition 

In the above calculations it was assumed that a a was anti-hermitian. A 
simple modification of the scalar product into a "rotated" one, allows to 
conclude that the equations are still symmetric hyperbolic (w.r.t. the new 
scalar product) even when the soldering forms are not anti-hermitian, as 
long as they can be transformed into ones that are so by a SL(2, C) trans- 
formation, that is, as long as the metric they generate is real and positive 
definite. 

Indeed, given a a such that the metric it produces is real, there exist 
SL(2,C) transformations L{a a ) such that a a := L~ 1 o a L is anti-hermitian. 
(See Appendix [] for a procedure to construct one.) Using any one of these 
transformations we define the following scalar product: 

(«2,«i)h ( «) = (H-WfaA^HWfaA 1 )) (9) 
= tr((L- 1 ^L)t L- x a\L) q ab + tr^IT 1 ^ L' 1 AIL) q a \ 

(10) 

where the linear operator H(u) is defined through L as in the second step. 
Thus, when computing 

P12H + P ] i2H = tr((L- 1 u 2 L)t J L- 1 PniL) + tr({L- 1 Pu 2 L)^ L^mL), 

we obtain terms with the structure,^] 

tT({L- 1 u 2 L) ] L- 1 {LaL- 1 )Ltr(u l LaL- 1 )) = 
tr((L _1 u 2 ^) t o-)tr(L _1 niL(7), 

that is, with the same structure as the one in the calculation in ||, and 
above, if one substitutes u%, and u 2 by their rotated versions, L~ 1 uiL, and 
L~ 1 u 2 L, and correspondingly a c by a c . Since a c is by construction anti- 
hermitian, symmetric-hyperbolicity follows for the extended system from 
the same calculations as in the previous subsection. 

As argued below, symmetric hyperbolicity in the reality condition sub- 
manifold does not seems to be enough for the usual proof of well posedness 
to work. Thus, a further modification of the equations is needed to obtain 
a set of equations which are symmetric hyperbolic in a whole neighborhood 
of the reality conditions sub-manifold. We do this in the next section. 



5 We consider the case N a = 0, for the part proportional to it is diagonal, hence sym- 
metric, without any conditions on a a . 
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3 Flow Regular izat ion 



Although, as we shall prove in section [5|, the evolution respects the reality 
conditions (and therefore solutions whose initial data sets satisfy them, have 
a metric which stays real along evolution), this property is not sufficient to 
conclude that the system is well posed. Indeed, all schemes used to prove ex- 
istence of solutions rely upon sequences of intermediate equations and their 
corresponding solutions, which are shown to form a contractive map, and 
hence to converge to the exact solution. These intermediate equations are 
in general linear versions of the original equations, but with the principal 
part evaluated at earlier approximated solutions. The problem is that this 
intermediate solutions do not propagate correctly the reality conditions and 
so one has to solve approximated systems in a neighborhood of the reality 
condition sub-manifold, but there the system is not symmetric hyperbolic, 
and so its eigenvalues have non-zero real part. Thus, the intermediate equa- 
tions unstable, and so if intermediate solutions exist at all their norm can 
not be properly estimated. 

This problem is overcome by further modifying the evolution equations 
outside the reality condition in such a way to obtain a symmetric hyper- 
bolic system even outside of it. Then standard contracting maps schemes 
for proving existence of solutions can then be applied to obtain local well 
posedness. 

Here we present a detailed description of the required extension to the 
flow in order to obtain a symmetric hyperbolicity in a whole neighborhood 
of the reality condition sub-manifold. 

If the metric is real and positive definite, that is, if it lies in the real- 
ity condition sub-manifold, then to prove symmetric-hyperbolicity one can 
simply use the fact that there is an SL(2, C) transformation which makes a 
anti-hermitian. Outside the reality sub- manifold no extension of the trans- 
formation would give an anti-hermitian soldering form, so the system there 
is not symmetric-hyperbolic without a further modification. 

The extra modification of the evolution equations consists first in extend- 
ing the SL(2,C) transformation outside the reality condition sub-manifold, 
with properties which ensure that when the reality conditions are satis- 
fied the resulting a a would be anti-hermitian; we present in detail one way 
to achieve this in Appendix |6| (the procedure is not unique). The second 
step is to change everywhere in the principal symbol of the system a a by 
La a L~ 1 := \(L{L~ l a a L - (L- 1 a a L)t)L~ 1 ). At the reality sub-manifold we 
re-obtain a a but, outside of it, the equations differ. Notice also that a a is, 
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by construction, anti-hermitian. Therefore we claim that the system is ev- 
erywhere symmetrizable, with symmetrizer H(a) := (L _1 )^L _1 to the left 
and H~ l to the right (notice that H is symmetric and positive definite ev- 
erywhere, since L £ SL(2, C)), that is, with respect to the following scalar 
product: 



<«2,«1>H (U ) = (i^- 1 (a)(a 2 ,^ 2 ), J ff(a)(a 1 ,A 1 )) 

= tr((L- 1 a 2 a L)t L~ l ~a\L) q ab + tr {{L~ 1 A 2 a L)^ IT 1 A\h) q ab , 

where q ab is the metric constructed out of a a , and therefore real. We shall 
denote with a ( " ) tensors which are constructed with a. Since they are 
SL(2,C) scalars they are also the ones constructed with Lo a L~ x . On the 
extended equations we shall raise and lower indices with q ab and its inverse. 

Again, we consider just the case N a = 0, for the part proportional to it 
is diagonal, hence symmetric, without any conditions on a a . Recalling the 
expressions for the principal symbols, 



P CT (a4ka)oi b = --mr{a 1 b a e )s a ec a c , (11) 
a 

P A (a,ik c )A\ = ^N^tr^ae)^ 6 ^ 

+ ^ e cdeh db tr (a e A d ) + -N e b dc a c tr(a a k [d Al)(12) 

Once modified, they become 

P <T (v,ik a )* 1 b = -^m?{d\La e L- l )e a ec La c L-\ (13) 

and 

P A (a,ik c )A\ = ^m [a tr(A 1 b] La e L~ 1 )i ae c La c L- 1 

+Ne cde ^La b L- 1 tr^a^ -1 A l d ) + e b dc Lo c L- x 
a a 

xtviLa-L' 1 k [d A a] ). (14) 

Which is identical to the one used in the previous section, but now a c is by 
construction anti-hermitian even outside the reality condition sub-manifold, 
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therefore La a L can be trivially rotated into an anti-hermitian one, and 
so symmetric-hyperbolicity follows as before. With this modification the 
standard proof of well posedness of symmetric-hyperbolic systems apply. 

Main Result: If we assume that data, (aft, A®), is given at a initial surface 
So such that it belongs (locally) to the Sobolev space H s (T i q) x H s ~ 1 (Eo), 
s > 3, and such that q ab and its first time derivative at Sq are very close to 
be real, ( that is, we start with initial data in a sufficiently small neighborhood 
of the reality sub-manifold) then a local solution exists and it stays in these 
spaces along the generated foliation. 



4 Characteristic Directions and Strong Hyperbol- 
icity 

It is of interest to see whether there are simpler modifications to the equation 
system outside the constraint sub-manifold which would give us well posed- 
ness (strong-hyperbolicity) but with a system not necessarily symmetric- 
hyperbolic. In doing this we shall calculate the complete set of eigenvectors 
with their respective eigenvalues. 



Recalling the expressions for the principal symbols of the system, (11), 
(|j~2|), we see that it is convenient to use at each point an orthogonal triad 
{k a , m+, m~} since we want to extend the results to a neighborhood of 
the reality condition sub-manifold. Let k a -the wave vector- be real, with 
|^a^ a | = 1 and let the other two vectors are taken to be null and orthogonal 
to k a . They are normalized such that, m^m^q ab := A;" 1 := (\/k a k a )~ 1 , 
and chosen so that when the metric becomes real m„ = (complex null 
vectors). Note in particular that we have the relation: 

£ a bc = i6&[ a m^m~j 

Thus the symbol can be diagonalized in blocks with the following eigen- 
values 

\ = ik a N a X + = i(k a N a +Nkq) A_ =i(k a N a -Nkq) 

and the subspaces associated with the above eigenvalues are 

E a ° = Span{fc a k d q d ,m+ k d q d ,m~ k d q d }, 
E a + = Span{/c a m~ d q d , ra+ m~ d q d , m~ m~ d q d }, 
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E c ~ = Span{k a m+q d ,m+ m+q d ,m- m+q d }, 

E A ° = Span{fc a k d q d , k a m+q d , k a m d q d }, 

E A + = Span{m+ k d q d ,m+ m+q d ,m+ m d q d }, 

E A - = Span{m- k d q d ,m^ m+q d ,m- m d q d }, 

In order to clarify how we have computed these eigenvectors and eigenval- 
ues, let us check this result for A\ = m^l d k d € E\, with l a G {k a , ?rt+, m~}. 
Consider the case N a = 0, for the part proportional to it is diagonal, then 

P A (a,ik c ) A\ = ?-m [a tT{Al]°e)e ae c d c 

+Ne cde ^a b tr(<r e A d ) + -N e b dc a c tr(a a k [d A a] ) 

= ctN (2k [b m+l e e ae c a c + 2k [a m+l d e a bc a c - k a m+l e e ade a b ) 

= -^aN (2e ban m +n e ae c l e a c + 2e adn m +n e a bc l d a c 
-e adn m +n e ad %a b ) 

= ifccrN mp d k d . (15) 

Note that, assuming the metric to be real and choosing the space- like 
foliation such that k a becomes the normal to a time-like hypersurface and the 
shift vector tangent to it, then we see that for each a a and A a , there are three 
incoming characteristics, the same number of outgoing (as should be the 
case for gauges respecting time-direction symmetry) and three characteristics 
which move along the boundary. 

Note that the eigenvectors span the whole space on the solution mani- 
fold, furthermore they are smooth functions on all arguments for k a ^ 0. 
Furthermore, if the metric is real, that is, if we are in the reality con- 
dition sub-manifold, then the eigenvalues are purely imaginary, thus we 
have a strongly-hyperbolic system. Can one make a simpler flow regular- 
ization and get strongly hyperbolicity (a weaker condition than symmetric- 
hyperbolicity)? That is, can we get away with the L transformation? From 
the form of the eigenvalues it is clear that this is the case. Indeed, one can 
modify the system outside the reality condition sub-manifold in such a way 
that the eigenvalues are purely imaginary also there, but since in general 
k is complex, and depends on k a , the needed modification transforms the 
differential equation system into a pseudo-differential one. We do not pur- 
sue this further because such a modification would not be practical for most 
applications one envisions. 
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5 The evolution of constraints and reality condi- 
tions 

Once a suitable extended evolution system for general relativity is shown 
to be symmetric hyperbolic in a whole neighborhood of the constraint and 
reality conditions sub-manifold, then standard results can be applied and 
different results on well posedness follow, in particular local evolution is 
granted if the initial data is smooth enough. For instance, for the specific 
system at hand, if the initial data (a a , A a ) is in .fP(E ) x H S ~ 1 (E ), s > 3, 
then the solution remains for a finite time in the corresponding spaces along 
the generated foliation E(. f] 

Given a solution in a foliation, we can identify via the lapse-shift pair 
any point on Ef for any t with a point at the initial surface, namely the 
points that lying in the integral curve of the four dimensional vector field, 
t a = Nn a +N a and so we can pull back to the initial surface the pair (a a , A a ), 
thus in the initial surface a solution can be seen as a one parameter family 
of fields in £P(E ) x fF-^Eo). 

The following geometrical picture of evolution emerges: we have an infi- 
nite dimensional manifold, K, of pairs of soldering forms and connections in 
a three dimensional manifold, Eo belonging to H s (T,q) x H s ~ 1 (T,q). Given 
any lapse-shift pair we can generate an evolution, that after the pull back is 
just a one parameter family of pairs in K, that is, a curve on that manifold. 
The tangent vector to that integral curve is our twice modified evolution 
equation system. Of course not all these integral curves are solutions to 
Einstein equations, for they would not generally satisfy the constraint nor 
the reality conditions, which in fact form a sub-manifold of K, denoted P. 

The relevant question is whether the integral curves which start at P 
remain on P, for they conform the true solutions to Einstein's equations. 
This would happens if and only if the tangent vector fields to the integral 
curves are themselves tangent to the sub-manifold P. In order to see that 
one can proceed as follows: 

We let RE, and CE denote the reality conditions quantities, namely 
^sq ab , ^STT a b, and the constraint equations quantities, namely C, C, and C a 
respectively. \\ We smear these expressions out with smooth tensors and so 



6 For the local problem one takes (through the choice of lapse-shift pairs) the foliation 
to be such that all the constant time surfaces coincide at their boundaries, that is, a Jens 
shaped domain of evolution. 

7 This is basically a coordinatization of a neighborhood of P in K. 
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obtain maps from the manifold K into the complex numbers, 

RE f Ja,A) := J f ab Zq ab dx 3 (16) 

RE fab (a,A) := J^f ab 5ir ab dx 3 (17) 

CE f -(a,A) := -i-v/2 f tr(fC) dx 3 (18) 

CE f (a,A) := -iV^ffCdx 3 (19) 

CE fa (a t A) := -iV2j^f a C a dx 3 , (20) 

We thus see that the reality-conditions-constraint sub-manifold P can be 
defined as the intersection of zero level set of each of this infinite set of maps. 
If these maps are sufficiently differentiable so that their gradients are well 
defined, it is clear that the tangency condition of the evolution vector field is 
just the requirement that when contracted with the differentials of the above 
defined maps the result should vanish at points of P. But these contractions 
are just the smeared out version of the time derivatives of the reality and 
constraint equations. Thus, provided we have enough differentiability, 8 we 
only have to check that this time derivatives vanish at P. 

Since these calculations need only be done at points of P, it only involves 
the original evolution vectors field, and so the standard results can be used, 
namely that the time derivative of the constraint equations, and the time 
derivative of the reality conditions vanish at P. The first result (see ap- 
pendix^) follows from the constraint algebra calculated in |1C]. The second 
one, basically, follows from the fact that the Ashtekar system is equivalent, 
up to terms proportional to the constraints, to Einstein's evolution equa- 
tions. Since in that equivalence one does not use any hermiticity nor reality 
condition explicitly, Ashtekar's equations are in fact equivalent to complexi- 
fied gravity that is to equations identical to Einstein's, but where the metric 
and the second fundamental form on each slice can be complex. Thus, at 
points where both the metric and the second fundamental form are real, 
and the constraints are satisfied, that is, at P, and provided that the three 
metric is invertible, the imaginary part of these tensors clearly vanishes. We 
have also verified this directly (see appendix ^, and see also previous works, 
Q and U), obtaining 

C t lTab = £ ^ab + N-KTT ab - 2Nir a e Tr eb - NR ab - D a D b N 



8 Note that in this manipulations no surface term arises, for in treating the local problem 
the lapse-shift pair vanishes at the boundary 
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-yft* (R + tt 2 - TTdcTT cd ) ~ iNe ab d D c (vr dc - Trq dc ) 
-2N ( -K d a K[db\ ~ nir[ ab ] + TT[ a d]Tr d b) - iNe b m D a 7r[ dm ] 



(D c Nir [md] -ND c ir [md] ) (21) 
where 

7T a fc = -tr(vr a a b ) and 7T[ a6] = tr(<T c P d cr d ). 

Since the terms, which could give an imaginary contribution to the expression 
when q ab and 7r a fc are real, are proportional to the constraints, we see that 
at points of P the evolution of the imaginary parts q ab , and ir ab vanishes. 
Note that the real part of ir ab corresponds to the extrinsic curvature, i.e. 
7T( afe ) = K ab . We conclude that, 

Main Result: Initial data satisfying both the reality conditions, and the 
constraint equations have an evolution which stays inside P , and so are 
solutions to the real and complete set of Einstein's equations. f\ 

We will conclude with a brief discussion of some issues relevant to nu- 
merical relativity. 

What does this mean when one does not solve Einstein's equations ex- 
actly, but just approximate them via numerical simulations or other means? 
If the approximations were a contraction map as the one often used to prove 
existence of exact solutions of symmetric hyperbolic systems, then it would 
follow that the approximate solution must approach the manifold P as it 
is refined and made closer to the exact solution. In practice the approxi- 
mation schemes are not contractive as required. Even after refinement, if 
the method yields convergence to the manifold P, this convergence could be 
very slow. 

To explore this problem, and in analogy with a similar work of Fritelli 
[13], we have looked at the evolution equations that the constraints satisfy 
when the flow is extended , but assuming that the reality conditions hold. 
From this study we have found that the constraint quantities satisfy by them- 
selves a symmetric hyperbolic system of equations. Thus, initial data sets 
which satisfy exactly the reality conditions, but are just near the constraint 
sub-manifold, if evolved by a scheme that respects the reality conditions 
then would stay "near" the constraint sub-manifold, in the sense that their 
departure would be bounded by a constant that depends only on the Sobolev 



9 The above argument is only valid for establishing local well posedness. For initial- 
boundary-value problems, one must be aware of boundary terms and impose conditions 
for them to also vanish. 
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norm of the initial data set. Unfortunately symmetric hyperbolicity by itself 
does not allow a finer control on the bound, which in principle could, even 
for linear equations, grow exponentially with the evolution time. Thus, this 
crude bound is not enough for controlling numerical simulations, but its 
absence would certainly make simulations very hard to implement. 

The evolution equations for the constraints have been calculated in ap- 
pendix ||. In order to prove the symmetric hyperbolicity we just need to 
consider the principal part of this set 

t = N a d a C-^=N[a a ,d a C], 

C = N a d a C - iV2Na 2 d a C a , (22) 

C a = N b d b C a + ^Nl a bc d b C c + ^=Nd a C. 

a 2 V2 

Then we prove the following, 

Lemma III.l: The equation system for any fixed, but arbitrary lapse 
and shift fields is a symmetric hyperbolic system in the reality condition sub- 
manifold. 

Proof: 

The eigenvectors and the eigenvalues can now be easily calculated, note 
that the principal symbol in this system 

P(y,ik a ) = P c ,(^,ik a )eP d (^,ik a ), (23) 

herey = (q,A), and C = (C, C b ). 

As above we use the orthogonal triad {k a , m a ,fn a } being m a and m a null 
vectors. Thus the symbol can be diagonalized in blocks with the following 
eigenvalues 

A = i k a N a A_ = i( k a N a - Nq) A+ = i( k a N a + Nq), 

and the subspaces associated with the above eigenvalues are 

E d ° = Span{/c d ^}, Eq~ = Span{ m d q d }, E d + = Span{m^ d }, 

E c + = Span{(0,m a ), (q, -^ k a)}, 

Eq~ = Span{(0,m a ), (q, --j=k a )}. 
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This set is a complete orthonormal set of eigenvectors with respect to 
the inner product 

(u^u 1 ) = {{C\C 2 ),{C\C 1 )) 

= tr(& 2 C 1 ) + -^C&C 1 + q ab C%}Cl 

and since the eigenvalues are purely imaginary the principal symbol is anti- 
hermitian. This concludes the proof of the Lemma.l(k 



6 Conclusions 

We have studied several aspects of Einstein's evolution equations as given in 
Ashtekar's formalism. We first have shown that when the reality conditions 
are satisfied, that is, when there is a SL(2, C) transformation that rotates the 
soldering form into a anti-hermitian one, the system is symmetric hyperbolic. 

Contrary to what one might expect at first , this condition is not enough 
to conclude that there are solutions, even local ones, to the evolution equa- 
tions. The problem being that, since the reality conditions are not linear 
conditions on the variables on which the problem is formulated, the usual 
contractive map of successive approximations does not respect them, and so 
evolution occurs in this approximations along paths where the system is not 
symmetric-hyperbolic, meaning that they can not be appropriately bounded 
by the usual methods. 

To remedy this problem we have proposed two further modifications. 
One of the modifications, called the regularized flows is a standard trick 
commonly used to grant symmetric hyperbolicity, suitably extended to the 
case under consideration, that is, to account for the intermediate SL(2, C) 
rotation mentioned above. In that case one obtains a symmetric hyperbolic 
system of equations even outside the reality condition sub-manifold. This 
allows then to establish the local well posedness of the problem by standard 
procedures. 

The second modification has the advantage of keeping unmodified the 
eigenvalues-eigenvectors structure of the principal part, which have a simple, 
and therefore probably useful, expression in terms of the solution, at the 
expense of transforming the differential equations into non local ones, namely 
into pseudo-differential equations. This modification can be implemented in 
numerical schemes -using fast Fourier transform- but only for a limited type 
of boundary conditions. 
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After the regularization of the flow and the subsequent implication of 
the existence of solutions of given differentiability, we discuss the problem 
of making sure that the solution, whose local existence we are asserting, 
does satisfy the whole set of Einstein's equations. In other words, whether 
integral curves to the evolution vector fields starting at the reality-constraint 
sub-manifold stay there. Since we have already shown that the solutions are 
unique (if sufficiently differentiability is assumed) then it is enough to show 
tangency of the vector field with respect to that sub-manifold. Provided 
some differentiability conditions are satisfied, this is equivalent to standard 
results which were nevertheless reviewed. 

We concluded section || with some considerations about the problem 
of numerical simulations, or approximations in general. Since errors are 
unavoidable in numerical algorithms one can at most consider evolution 
in a hopefully small neighborhood of the reality-condition-constraint sub- 
manifold. Thus, even when -in the best of the cases- knowing that eventu- 
ally, by refinements of the approximation scheme, the solution would con- 
verge to that sub-manifold, one would like to have a priori estimates of the 
deviation as a function of the initial error. 

We have shown, following the ideas of Fritelli for the ADM formulation, 
that the constraint system -while at the reality sub-manifold- satisfies by 
itself a set of symmetric hyperbolic equations. Thus such an a priori bound 
follows directly. We consider this a preliminary result, for the bounds that 
follow from symmetric hyperbolicity are too crude for numerical purposes, 
and so a further refinement, this time using information about the full sys- 
tem of equations (contrary to symmetric-hyperbolicity which only uses the 
principal part structure of them), seems to be needed. 

We will conclude with a remark. A substantial improvement on the 
handling of these equations would be the possibility of imposing reality 
conditions in a linear way, that is, by requiring the hermiticity (or anti- 
hermiticity if considered as matrices) of the soldering forms. Unfortunately 
with the present scheme of modifying the equations outside the constraint 
sub-manifold, the hermiticity condition restricts unnecessarily the possible 
evolutions, for it basically fixes a unique lapse. Otherwise the system can not 
be made symmetric hyperbolic in an straightforward way. This restriction is 
unnatural, for one would expect that this condition could be enforced fixing 
the SU(2) gauge freedom, and not the one associated with diffeomorphisms. 
But all attempts to fix the SU(2) gauge appropriately seems to involve (at 
best) elliptic conditions on the "time component" of the connection. Is there 
an alternative avenue to handle the linear reality conditions? 
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APPENDIX A: THE CONSTRUCTION OF THE SL(2, C) 
TRANSFORMATION. 

Given a soldering form, a a , we shall define here a procedure to construct a 
SL(2, C) transformation such that, when the metric is real the transformed 
soldering form is anti-hermitian. Given a transformation L making a a anti- 
hermitian, there is a whole set of them with the same property, namely one 
can left-right multiply the new one by a SU(2) rotation and obtain a new one 
with the same anti-hermiticity property. To remove this arbitrariness in our 
construction we fix a real basis (not orthonormal), {ef} i = 1,2,3 of vector 
fields. For simplicity we shall work here with frames, {-Ef} i = 1, 2, 3 so that 
a a = Efr % where t % is an anti-hermitian Pauli set of matrices. We then have 
EfEjS^ = q ab , and so EfE^q^ = <5jj. We define L as the transformation 
such that the transformed frame {Ef} i = 1, 2, 3 satisfies: 

Ef = b\e\ 

E a 2 = b\e a l+ ble a 2 + blel with |, | real. 

The ortho-normality condition then fixes the third vector, and so the com- 
plete transformation. Notice that if the initial frame can be rotated to a real 
one, that is, if the metric is real, then the rotated frame we are constructing 
will be real. 

We let Ef = Cje% and Ef = LjE* = WCj k e a k = Cje a k . Since is 
invertible, knowing C is equivalent to knowing L, we shall construct C. 

We have, C\ l = b\, Ci 2 = £i 3 = 0, b\ is fixed by the ortho-normality 
condition: 1 = (b\) 2 ei ■ e±, we choose the root with positive real part. 

Next we have C2 1 = b\ C2 2 = b\ C2 3 = b^, and two conditions, the 
normality condition, 

1 = (bl) 2 e 1 -e 1 + {b 2 2 ) 2 e 2 -e2 + {bl) 2 e 3 -e 3 + 2blble 1 -e2 + 2b 2 2 ble2-e3 + 2blble 1 -e 3 . 



and the orthogonality one, 

= b\(b\ei ■ ei + b\e x ■ e 2 + b\e\ ■ e 3 ) 
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We require that 9(§jtff) 0) otherwise we impose the reality condition 
to 62/62 an< ^ take b\ = 0. We multiply the orthogonality relation by a phase 
e such that $s(e i0 ei ■ e^) = 0, we then divide by b\ and assume the quotient 
of the fr's to be real. Taking real and imaginary part we get two real linear 
equations: 

|9(e*%- ei ) + |9( e *%-e 3 ) = 

§H(e* ei ■ ex) + §^(e l % • e 3 ) = -|ei • e 2 |, 

which, because of the requirement above, have a non-vanishing determinant. 
Thus, it has a real solution for the quotients. We then solve the normality 
condition for b\, choosing, as before, the root with the positive real part. 
The remaining coefficients of C are completely determined by the remaining 
ortho-normality conditions (three equations) and the condition that the real 
part of 63 be positive. 



APPENDIX B: BRIEF OVERVIEW OF THE HAMILTON! AN 
FORMULATION IN ASHTEKAR'S VARIABLES 



In this section we shall present a brief overview of the Hamiltonian formula- 
tion in spinorial variables in order to get Einstein field equations as a system 
of evolution and constraints equations and thereby to study the evolution of 
the constraints (we follow [10| and [111]]). 



1 The Lagrangian Framework 

Let us consider a four-manifold M which has topology SxR, for some three- 
manifold X, with a four dimensional SL(2,C) soldering form, o a a A > an d a 
connection 4: A a A B which acts on the unprimed spinor indices. We restrict 
the soldering form to be anti-Hermitian so that it defines a real space-time 
metric via g ab = (^ AA ,a bAA with signature ( — h ++) and a unique torsion- 
free derivative operator V compatible with a a A A ■ Finally, the derivative 
operator A V a defined by 4 A a A B via 4 P q Aa = <9 q Aa + A A a A B \B acts on 
unprimed spinor s. 

The gravitational part of the Lagrangian density of weight 1 is: 



r (A„\„a A'b 4p 

L = { a)a a & BA> ^ab 
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where ( 4 cr) is the determinant of the inverse soldering form and A F ab A B is 
the curvature tensor of 4 A a A B - 

It is useful to perform a 3 + 1 decomposition of the action, this shall be 
needed latter in order to pass on to the Hamiltonian framework. 

Let us introduce on M a smooth function t whose gradient is nowhere 
vanishing and whose level surfaces S t are each diffeomorphic to S. Let t a 
be a smooth vector field on M with affine parameter t and on each level 
surface, let n a be the future-directed, unit, time-like vector field orthogonal 
to S t . Denote the induced positive-definite metric on T, t by q ab = g ab + n a n b 
and obtain the lapse and shift fields N and N a by projecting t a into and 
orthogonal to E t , i.e. t a = Nn a + N a . 

Identifying unprimed SL(2, C) spinors on M with SU(2) on £ t , we intro- 
duce the soldering form a a a B on SU(2) spinors. It defines the three metric 
on X( as q ab = — tr (a a a b ). Here a matrix notation is employed, and shall be 
used in the following, for unprimed spinor indices in which adjacent summed 
indices go from upper left to lower right, e.g., (a a a b )A B = cr a A C o'c B - 

Let A a A B , F a bA B and V a (the Sen connection) be the pull-backs to S t 
of 4 A aA B , A F abA B and 4 V b respectively. Then V a \ A = d a X A + A aA B \B and 
F ab = 2d [a A b] + [A a ,A b }. 

Finally, there is a natural (canonical) spinorial connection associated 
with the three-metric such that D a a b = 0. It relates to the Sen connection 
via the extrinsic curvature K ab on E t as follows 

F> a ^A = E> a \A + -j=K a A B \B 

= d a \A + ^aA B ^B + -^K a A B \B, 

where K a A B = K ab a a A B and T a A B is the spin connection 1-form of D. Then 
A a = T a + -^K a and using the fact that the derivative D is compatible with 
a a , i.e. 

D a o- b = d a a b - T ab c a c + [r o , a b ) = 
with T ab c denoting the Christoffel symbols; we calculate 

e bcd a a e bcd a b 
Ta = ~2Vf tr t^ 9 ^ VT tT ( aadc(Td ^ ^ 

where the orientation three-form on T, t is written as e abc = — y/2ti (a a a b a c ). 
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The action Q can be expressed in terms of only of three-dimensional 
fields: 



S = J dt J d 3 xti(V2ia b (£ t A b -V b ( 4 A-t))-N a a d b F ab -V2iN a d h F ab 

where a a ab = (o") °~ a AB, N = (<t) _1 N and the Lie derivatives treat internal 
indices as scalars. 

The action depends on five variables N, N a , 4 A ■ t, A a A B and a a a B ■ 

The first three variables play the role of the Lagrange multipliers, only the 
last two are dynamical variables. Varying the action with respect to the 
Lagrange multipliers we obtain the constraint equations: 



C(a,A) 
C a (a,A) 
C A B (v,A) 



ti{d a a b F ab ) = 0, 

ti{a b F ab ) = 0, (25) 
V a a a A B = 0, 



and varying with respect to the dynamical variables, yields the evolution 
equations: 

C t a b = -l 4 A.t,a b ] + 2V a (N^) - ^=V a (N[a a ,a b }), 

C t A b = V b ( 4 A-t) + N a F ab + ^=N[a a ,F ab }. (26) 



2 The Hamiltonian Framework 

Recall that the dynamics of a mechanical system can be achieved having 
what is called a symplectic manifold, i.e a pair (r, O a fl), where T is an even- 
dimensional manifold, and Q a p a symplectic form, i.e a 2-form which is closed 
and nondegenerate. Given any function / : V — ► R, the Hamiltonian vector 
field of / is defined by Xf a = £}P a Vpf. Given any vector field v a S T X T, 
we say that v a is a symmetry of the symplectic manifold if it leaves the 
symplectic form invariant, i.e if C v VL a p = 0, in which case the diffeomorfisms 
generated by v a are called canonical transformations. 

Given two functions /, g : T — * R, their Poisson bracket is defined by 

{f,g} := Q al3 V a fVpg = C Xf g = -C Xg f. 



The surface terms are not included here. 
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Thus the dynamics of a physical system is given by assigning a phase space 
(r, £l a s, H) on which the evolution is generated by the Hamiltonian H from 
the initial state. Hence for any observable / 

f = C X J = {H,f}. 

When the system is constrained, i.e. when there are points in the phase 
space that can not be reached by the physical system, the system remains in 
a sub-manifold called constraint sub-manifold, that can be specified by the 
vanishing of a set of functions T = {p £ V /Ci{p) = 0, for i = 1 • • • m}. The 
constrained system is said to be first class if there exist functions , i,j,k = 
1, • • • , m, called structure functions, such that {Cj, Cj} = fij k Ck- 



In the formulation of General relativity in Ashtekar variables, the config- 
uration space C is the space of all weighted soldering forms a a , and the phase 
space T is the cotangent bundle over C. The phase space is represented by 
the pairs (a a ,A b ) (these variables, apart from a numerical factor, the new 
canonically conjugate pair, ). The action of a cotangent vector A a on any 
tangent vector [5a) a at a point a a of C is given by 



A(S<r) = - J d 3 x tr(A a 5a a ). 



We will not discuss boundary conditions here, but just remark that one 
possible choice would be to require that the canonically conjugate fields 
should admit a smooth extension to the point at spatial infinity if the three- 
surface £ is made into the three sphere by the one point compactification, 
for an extensive discussion of fall-off properties of the fields see flofl . 

Thus in order to obtain the constraint sub-manifold where the physical 
gravitational states take place, we need to construct functionals from the 
constraints (^), i.e. we need to smear out these constraints with a function 
N, a vector field N a and an anti-hermitian traceless iVy^ test fields. We 

define 

C N (a,A) := -iV2[d 3 x tr(NV a a a ), 

C N {a,A) := -iy/2 f d 3 xN ti{d a d b F ab ), (27) 

C^{a,A) := -iV2 f d 3 x N a tr{a b F ab - A a V b a b ), 

then we should ask if these constraint functions are first class in the above 
terminology, to answer this question we have to calculate their Poisson brack- 
ets. For that purpose, we introduce the symplectic form 

n\(a,A) ((6a t 5A), (S^M')) ■= ^= / s d 3 x tx{8A a 5a' a - 5A' a 5a a ), 
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where (Sa, 5 A) and (5a', 5 A') are any two tangent vectors at the point (a, A). 
Thus given two functionals / and g their Poisson bracket is 

r , , * f ,3 ( Sf Sg Sf Sg \ 
{f > 9} = T2h dXtr {^^-^SAa)- 

The Poisson bracket between the constraint functionals and our fundamental 
variables are: 

{a a MN ,A b AB } = -±=5(x,y)8 M {A 8N B \ 

{C N ,a a } = [a a ,N], 

{C N ,A a } = V a N, 

{C N ,a a } = 2V b (Na^ a a\ (28) 

{C N ,A a } = N[a b ,F ba ], 

{ C N^ a } = -^n^, 
{Cf[,A a } = —C^A a , 

and the constraints algebra becomes 



{Cn,Cm} = C[JV,M]> 



{Cfi,C M } = c Cfi 



{^JV'^MJ — C[JV,M]' 

{C N ,C M } = 0, (29) 
{ C N' C M} = C C S M 

where K a = 2a 2 (MV a N -NV a M) . 

Using the canonical variables (a a ,A a ) and the action, the resulting 
Hamiltonian writes: 

H = d 3 x ( i V2tr(N a a b F ab - N a A a V b a b ) +N tr( a a a b F ab )^j , 

= -C$ + ^C N , (30) 

where the Lagrange multiplier A A a t a has been chosen as A a N a . 
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Therefore, the Hamiltonian evolution of the dynamical variables, yields: 
Cfa b = {H,a b } 



-[A a N a ,a b ] + 2V a (NK b ^+N b V a d a - -^=V a (N[a a , a 6 ]) (31) 



L t A b = {H,A b } 



C^A b + ^=N[a a ,F ab ] 



= V b (A a N a )+N a F ab + ^=N[d a ,F ab }. (32) 

Note that the evolution obtained from the Lagrangian differs from the evo- 
lution obtained from the Hamiltonian in a " constraint term" . 



APPENDIX C: THE CONSTRAINTS EVOLUTION 



To compute the evolution equations for the constraints with the ex- 
tended flow let us first calculate the Hamiltonian evolution of the constraints 
(Ca B ,C,C a ). Using the Hamiltonian given by (p0|) , the constraint algebra 
( p9| ) and integrating by parts, we obtain 

t = {H,C} = C M C, 

C = {H,C} = C^C + iV2 a 2 N V a C a + i 2 V2a 2 C a V a N. (33) 
C a = {H,C a } = C^C a -^NV a C-iV2CV a N+^=Ntr([a b ,F ba }C). 



The most complicated calculation is the last one (the evolution of C a ) 
in order to do this, we define the functional 



= -n/2~y d 3 x M a ti(a b F ab ), 



then 

i 



{H^C^} = -{Cjf,C a } + -j={qf,C a }. (34) 



V2 

We calculate each one of the terms above as follows 



> 



{CfrCtf} = -= d d xtr 



6C N 6C M SC N SC M 



y/2 7e 1 SA a Sa a 5a a 5A a 
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= {Cjf,Cjtf} + {Cfl,Cp}\p=M«A a 
+iV2 [ d 3 x C$A a M a C, 

J S 

= C [N,M] + C C f} (M^A a ) + ^ C M «C s A a 

= C [N,M] ' ( 35 ) 

and 

{ C Ni C m) = { C Ni C m) + { C N C P}\P=M"A a 

-iy/2 [ d 3 xNtv([a b ,F ab ]M a C), 
= ~ C C lCl N+ C N[ab,F ab ]M a - ( 36 ) 



Inserting the Poisson brackets ( j35|) and (|3q ) in (34), using that the Lie 
derivative of a function with weight minus one is C^N = M a d a N — N d a M a 
and integrating by parts we get 

{H, C a } = C^C a - -^NV a C - iV2CV a N+ -l=NtT([a\ F ba ]C). 

Since we have changed the evolution of the dynamical variables outside the 
constraint sub-manifold, the calculus of the time derivative of the constraint 
equations must be done using equations (||) and (||). If / : V — > C then 

and since the perturbation we have made on the equations (||) and (B) are 
linear in the constraints, we obtain 



/(*, A) = {H,1}+ J^d 3 x 

Then, using (|33|) in (|37|), we have the constraints evolution 

C = N a d a C-^=N[d a ,V a C}-^=V a N[a a ,C]-i2V2NC d a d , 
C = N a D a C -iV2Na 2 D a C a -iV2C d Ntr(Ca d ^) 
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+-LNtr([C,a a }[a b ,F ab }) + C~e abe tv ((V b a a ) a e ) , (38) 
C a = N b D b C a + ^Ni a bd D b C d + ^=ND a C + C b D a N b 
--^C d D b N~e a db - ^Nti ((V [a a lcl ) a b ) ~e b \ dc C d . 

APPENDIX D: THE REALITY CONDITIONS 

Consider the variables (a a ,A a ). We want to prove that if they define 
a real metric on an initial surface, then this metric remains real under the 
Hamiltonian evolution. In order to do this, we shall calculate the evolution 
of the variables q ab = — tr (a a a b ) and T[ ab = — tr (ir a a b ), where Tr a is defined 
by 

D a \A = D a \A + —j=ir a A B ^B- 

We write the evolution of a a in terms of ir ab as follows 

C t a b = C^ b -^V a ^a^a b ^+^Np a a a ,a b ] 

= C n a b -i V C N a e cbe a e + aN (n b e a e - vra 6 ) , (39) 
then we calculate 

C t q ab = -2tr (Ctd^a^ 

= C^q ab + 2Na 3 (vr^ - nq ab ) . (40) 

The evolution of q ab follows from the evolution of a 2 = det(q ab ), thus we 
calculate 

C t a 2 = 2a 2 D a N a - 2a 3 irN. (41) 

Finally, we get 

£ t q ab = £fiq ab -2NTT( ab y (42) 

Thus in order to ensure the reality of q ab we need to know the reality of n ab . 
The evolution of Tr ab is given by 

C-tKab = -tr (n a £ t (Tb) + iV2tr (a b C t A a ) - iv^tr (<r b C t T a ) . 

Hence in order to calculate it, we need to rewrite the evolution of A b and a a 
in terms of ir ab . Using 

Fab = Rab ~ ^[a^b] + iV2D [a lT b] 
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where R a b is the spinorial curvature, and redefining C = tr ( o~ a a b F a b ) and 
Cb = tr (cr a F a b), we compute 



C t a h = C$a b -iD c N a e c be a e -Nir eb a e (43) 
C t A b = £^A b + -^=N[a a ,F ab ]--^=Na b C - iNe b dc a c C d 



V2 L ' a ° J 



i 



£,A b + -^=N[a a ,R ab } 



--j=N (vr an vr b a - vrvr^) a" - i\r[<7°, ^ [a vr fe] ] 
-^=Na b C-iNe b dc a c C d . (44) 

Hence 

£ t vr afe = C M ir ab + iV^vr^ - 2Nir eb Tr a e - NR ab - q ab NC + NV2e ab d C d 
-i2Ne b de tr (D [d ir a] a e ) + iD c N e ce b 7r ae 

+iV2ti ( (£ ^r a - C t T a ) o b ) . (45) 

Using the fact that D a a b = and the formula for r a , the last three terms 
can be set as 

-i2Ne b de ti (D [d TV a] a^ + iD c Ne ce b vr ae + i^tr ((c$Y a - C t T a ) a b 



-D a D b N - iNe b dm D a ir dm - -q ab e cmd (D c Nn md + ND c n md ) 



yielding 



CtKab = £ ff^ab + N7TTT ab - 2A% a e 7T e6 - N R ab - D a D b N 

'r + tt 2 - TT dc 7r cd ) - iNe ab d D c (7T dc - TTq dc ) + 



2 

2N (Tr d a ir [db ] - iriT [a b] + Tr[ ad ]ir d b 



-iNe b dm D a 7r [dm] - -q ab e cmd [D c N7r [md] - ND c n [md] ) . (46) 
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